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Preface 


This  report  is  a  compilation  of  presentations  given  at  the  ‘•Workshop  on  'I^ajcc- 
tory  Optimization  Methods  and  Applications”,  held  at  the  1992  AIAA  Atmospheric 
Flight  Mechanics  Conference  in  Hilton  Head,  South  Carolina.  This  workshop  was 
co-chaired  by  Harry  Karasopoulos  and  Kevin  J.  Langan,  both  of  the  former  Flight 
Performance  Group  of  the  High  Speed  Aero  Performance  Branch  in  Wright  Ljvbora- 
tory. 

It  is  hoped  that  this  document  will  help  the  attendees  retain  some  of  the  itlcas 
presented  in  the  workshop,  in  addition  to  providing  useful  information  to  those  who 
were  unable  to  attend.  Appreciation  is  expressed  to  the  presenters  and  attendees. 
For  the  third  year  in  a  row,  this  workshop  has  proved  to  be  a  successfxil  forum 
for  highlighting  current  work  in  trajectory  c^timization.  Thanks  are  also  dtic  to 
the  American  Institute  of  Aeronautics  and  Astronautics  for  making  this  workshop 
possible. 

The  following  was  the  workshop  sdbedule: 

SCHEDULE 


•  1:00  -  Introduction  -  Harry  Karasopoulos,  Wright  Laboratory. 

•  1:00  -  OMAT:  An  Autonomous  Optimal  Solution  to  Rendezvous  Problems 
with  Operational  Constraints  -  Don  Jezewski,  McDonnell  Douglas  Si)acc  Systems 
Company,  Houston,  TX. 

•  1:15  -  MULIMP:  Multi-Impulse  Trajectory  and  Mass  Optimization  Pro¬ 
gram  -  Darla  German,  Science  Applications  International  Corporation,  Scliamiil)iirg, 
IL. 

•  1:30  -  Phillips  Laboratory  Applications  of  POST  -  Jim  Eckmann,  SPARTA 
Inc.,  Edwards  AFB,  CA. 

•  1:45  -  OTIS  Advances  at  the  Boeing  Company  -  Steve  Paris,  The  Boeing 
Company,  Seattle,  WA. 

e  2:00  •  OTIS  Activities  at  McDonnell  Douglas  Space  Systems  Company  • 

Rocky  Nelson,  McDonnell  Douglas  Space  Systems  Company,  Huntington  Beach,  CA. 

Hi 


•  2:15  -  Advances  in  Trajectory  Optimization  Using  Collocation  and  Non¬ 
linear  Programming  -  Dr.  Bruce  A.  Conway,  University  of  Illinois,  Urbana,  IL. 

•  2:45  -  BREAK 

•  3:00  -  Flight  Path  Optimization  of  Aerospace  Vehicles  Using  OTIS  -  Dr. 
Rajiv  S.  Chowdhry,  Lockheed  Engineering  and  Sciences  Company,  NASA  LaRC, 
Hampton,  VA. 

•  3:15  -  Trajectory  Optimization  of  Launch  Vehicles  at  LeRC:  Present  and 
Future  -  Dr.  Koorosh  Mirfakhraie,  ANALEX  Corp.,  NASA  Lewis  Research  Center, 
Cleveland,  OH. 

•  3:30  -  Collocation  Methods  in  Regular  Perturbation  Analysis  of  Optimal 
Control  Problems  -  Dr.  Anthony  J.  Calise,  Georgia  Institute  of  Technology, 
Atlanta,  GA. 

•  3:45  -  Automatic  Solutions  for  Take-Off  from  Aircraft  Carriers  -  Lloyd  II. 
Johnson,  AIR-53012D,  Naval  Air  Systems  Command,  Washington,  DC. 

•  4:00  -  Current  Pratt- Whitney  OTIS  Applications  -  Russ  Joyner^,  United 
Technologies,  Pratt- Whitney,  West  Palm  Beach,  FL. 

•  4:15  •  Airbreathing  Booster  Performance  Optimization  Using  Microcom¬ 
puters  -  Ron  Oglevie,  Irvine  Aerospace  Systems  Co.,  Fullerton,  CA. 

•  4:30  -  Scheduled  Session  End 

Unscheduled  Speakers 

•  Dr.  Klaus  Well,  University  of  Stuttgart.^ 

•  Dr.  Mark  L.  Psiaki,  Cornell  University. 


‘Cancelled 

^Presentation  copy  not  available  at  the  time  of  printing 
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Analytical 


Numerical 
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Pefipitlon  Qf  an  Optimal  Bendfizvous  >?.roblem 

□  Given: 

■  Chaser  and  Target  States:  (Sc.tc).  (Sy.tT) 

■  Attracting  Body 

■  Objective  Function,  i.e.,  Oelta-V,  Fuel,  Time,  etc. 

□  Subject  To: 

■  Force  Field 

•  Perturbations 

■  Terminal  &  Inflight  Constraints  and  Limits 

□  Define: 

•  Optimal  Sequence  of  Maneuvers  (impulses,  Finite  Bums) 

•  Number,  Location,  Magnitude  or  Duration,  and  Direction 

□  Such  That: 

•  Chaser  &  Target  Vehicles  Achieve  a  Relative  Configuration  at  Some 
Time 


V _ 
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Extremal 

Solution? 


Compute  I 

x,;;  I 

Iterate  NZSOL 


Variation 
of  Solution 


ja,vj.v6,vc 


X  -  Parameter  Vector 
S  -  State  Vector 
A  -  Costate  Vector 
J  -  Ob|ective  Function 
G  -  Active  Constraint  Vector 
jC  •  Solution  Cost  ^ 
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Forms  of  Constraints 

Q  Three  Types  of  Constraints 

■  Parameter  Constraints 

• a<  X<  b 

■  Linear  Constraints  (Constant  Jacobian  Matrix) 

•  L{S,t)  ^  0 

■  Non-Linear  Constraints 

•  NL(S.t)  ^  0 
Q  Bounds 

■  Constraints  are  Bounded  by  Upper  and  Lower  Bounds  (Bi^.Bu) 

■  Equality  Defined  by  By  «  Bi, 

■  Unbounded  Defined  by  Bl»  -  «»,  Bu=  + «» 


^ _ 

AfeOo/wMtf  Oou0te«  Spae^  Syatama  Company -Houston  DMsIon 


Paga6ot22 


- - - 

Constraints  are  Defined  bv  Five  Integers 

□  1  -  Constraint  Number 

□  J  -  What  the  the  Constraint  is  Referenced  to: 

•  An  Impulse 

•  Another  Constraint 

□  K  -  Reference  Impulse  or  Constraint  Number 

□  L  -  Condition  that  Triggers  or  Initiates  Constraint 

•  Time  from  GMT  or  Reference  Event  (Impulse  or  Constraint) 

■  Phase  Angie 

■  Lighting  Condition  between  Chaser  and  Target  Vefiicies 

■  Delta-Angular  Measurement  in  Chaser  Orbit 
"  Number  of  Revolutions 

■  Chaser  Vehicle’s  nth  Periapsis,  Apoapsis  Crossing 


V _ t 
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-  : 

Constraints  are  Defined  bv  Five  Integers  fCQncJ’d) 

□  N  -  Type  of  Constraint 
•  Periapsis,  Apoapsis 

■  Differential  Height  between  Chaser  Vehicle  &  Target  Orbit 

■  Phase  Angle 

■  Sleep  Cycle  or  Quiet  Time 

■  Chaser  Orbit  Coelliptic  with  Target  Orbit 

■  Chaser  Position  Vector  Relative  to  Taurget  LVLH  Frame 

■  Chaser  Velocity  Vector  Relative  to  Target  LVLH  Frame 

■  Bounded  Delta*V  in  Chaser  LVLH  Frame 

■  Inertial  Line-of-Sight  Angular  Rate 

■  Wedge  Angle 


Vs 


UeDoiwll  Dougtaa  Spacm  Syatama  Company  •  Houaton  Dhrtalon 


Paoa8ol22 


Demonstration  2: 
Typical  Shuttle  Rendezvous 


□  Shuttle  (Chaser) 

■  110  circular  altitude 

•  Inclination  =  28.5’ 

•  Longitude  of  ascending  node  =  101’ 

■  Argument  of  perigee  =  0’ 

■  True  anomaly  =  180*  q 

□  SSF  (Target) 

■  190  nmi  circular  altitude 

>  Inclination  =  28.5’ 

"  Longitude  of  ascending  node  =  100’ 
"  Argument  of  perigee  =  0* 

>  True  anomaly  =  0’ 

□  Limited  to  2  maneuvers 


(Not  to  Seal*) 
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Demonstration  2  (Cont’d) 

□  Results 

■  OMAT  handies  perigee  constraint  for  unperturbed  orbits 

■  Optimum  unconstrained  trajectory  required  461  ft/s, 
constrained  trajectory  required  603  ft/s 
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Perturbations 

□  Primer  Vector  Theory  (Presently)  Requires  Perturbations  to  have  Form 

•  9t(R.t) 

□  Largest  Geopotential  Perturbation  for  Earth  &  Mars  Is  J2,  Has  Form 

^J2 

•  Where  K  is  a  Unit  Vector  Normal  to  the  Equator  and  and 
are  Functions  of  the  Position  Vector. 

□  Presently  Incorporating  NxM  Geopotential  Model 

□  Need  to  Extend  Theory  to  Functions  9i(R,V,t) 

□  Need  to  Develop  Theory  for  Third-Body  Effects  (Libration  Point  Rendezvous) 
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Demonstration  4;  I 

Mars  MEV  Post-Ascent  Rendezvous  with  MTV  | 

a  MEV  (Chaser) 

■  117  X 135  nmi  altitude 

>  Inclination  r  164.264662* 

■  Longitude  of  ascending  node  =  194.39079* 

>  Argument  of  perigee  s  159.820571 

■  True  anomaly  s  0* 
a  MTV  (Target) 

■  135  X  18,294  nmi  altitude 

■  Inclination  =  164.2* 

>  Longitude  of  ascending  node  =  195* 

•  Argument  of  perigee  s  16* 

•  True  anomaly  s  0* 
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a  A  Definition  of  Lambert's  Problem: 

■  What  is  the  Initial  Velocity  Vector,V|,  that  Generates  a  Trajectory  that 
Passes  between  Two  Radii  Separated  by  a  Given  Angle  in  a 
Specified  Time  Interval  ? 

V|_ 
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□  Standard  Approach 

•  V|  Obtained  from  Classical  Lambert  Problem 

•  5V|  Obtained  from  Solution  to  Variational  Equations,  i.e., 


<I>(t,.t,) 


5V.  =  <^-;6R,  <I»(t..t,,= 

^21  ‘>22] 

Difficulties  with  this  Approach 

•  5Rp  Must  be  ‘Smair  (Linear  Approximation) 

•  Conditioned 

•  J2  Frequently  Must  be  Reduced  (Sub-Problem) 

•  Excessive  Number  of  Iterations  and  Integrations 
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Improved  Approach 

•  V|  Obtained  from  First-Order  Correction  to  the  Inverse-Square 
Problem  Resulting  from  the  J2  Perturbation 

-  Analytic  Solution 

-  Expressed  In  “Ideal  Referen<»  Frame" 

-  Regular,  No  Singularities 

-  Solution  Expressed  in  Terms  of  Elements  &  their 
Variations 

•  5V|  Obtained  from  Solution  to  Variational  Equations,  i.e., 

6V,.<t.-;25R, 

•  5Rp  =  0(10-3)  Smaller  than  Classical  Approach 

•  No  Requirement  for  J2  Reduction  (Sub-Problems) 

•  Solution  to  TPBVP  Requires  Only  a  Few  Iterations  and 
Integrations 
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a  Unperturbed  Problem 

■  Force  Field  (Inverse-Square),  i.e.. 

V  -  -n(R/r^) 

■  State  and  Variation  in  State  Obtained  from  Solution  of 

•  Kepler’s  Problem  (Goodyear,  Analytic) 

■  Boundary  Value  Problem  Satisfied  by  Solution  of 

•  Lambert’s  Problem  (Gooding.  Analytic) 

■  Constraints  and  Variation  of  Constraints  Evaluated  Using 

•  Keplerian  Elements 

■  Non-Linear  Programming  Algorithm.  NZSOL,  Solves  Constrained 
Optimal  Problem 

■  Solution  Obtained  in  seconds  on  Sun  Sparcstation  2 
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Outline  of  Optimal  Solution  Approach  (concrd^ 

Q  Perturbed  Problem 

■  Force  Field  (Perturbed  Inverse-Square) 

V  =  -^(R/r^)>SR 

■  State  and  Variation  in  State  Obtained  by  Numerical  Integration  of 
Differential  Equations 

•  Variable-Step  Runge-Kutta-Fehiberg  7/8  Algorithm 

•  Maximum  of  42  Linear  &  Non-Linear  Differential  Equations 

■  Boundary  Value  Problem  Satisfied  by  Solving 

•  Perturbed  Lambert  Problem 

■  Constraint  Targeting  Evaluated  Using  Non-Periodic  Elements 

•  Constraints  and  Variation  of  Constraints  Evaluated  on  Perturbed 
Trajectory 

■  NZSOL  Solves  Perturbed.  Constrained  Optimal  Problem 

■  Solutions  Obtained  (Presently)  in  Minutes 


UeOonnM  Dougin  Spac*  Company  •  Houtton  DMghn  Paaa20ol22 


- - - - - 

WbaLfclave  We  Done.  Where  Are  We  Going? 

Q  Present  Status  of  Development; 

■  Proof  of  Concept.  Using  Unperturbed  Solution  to  Solve  Perturbed 
Problem 

•  Verify  Solution  Approach  for  Handling  Perturbations  and  Constraints 
and  their  Conditions 

■  Developed  Solution  Approach  for  Perturbed  Lambert  Problem 

■  Illustrate  Initial  Capability  of  the  Algorithm.  (OMAT).  to  Efficiently 
Solve  Optimal  Rendezvous  Problems  with  Operational  Constraints 

□  Planned  Future  Development; 

■  Expand  Perturbation  and  Constraint  Models 

■  Develop  Multi-Rev.  Capability 

■  Develop  Finite-Thrust  Model 

■  Libration  Point  Rendezvous 

■  Solve  Advanced  Problems 
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Condudina  Remarks 

□  Integration  of  Theoretical  &  Operational  Aspects  Of  Optimal  Rendezvous 

□  Development  of  an  Autonomous  Optimal  Rendezvous  Solution 

□  Primer  Vector  Theory  Basis 

□  Approach  Based  on  Solution  to  Lambert’s  Problem  and  it’s  Extention  in  a 
Perturbed  Force  Field 

□  Premise  is  Made  that  Perturbed  Problem  is  an  e  away  from  Unperturbed 

□  Constraints  are  Adjoined  to  Objective  Function  by  Lagrange  Multipliers 
Q  General  Mapping  for  Constraints  from  State  to  Parameter  Space 

□  NZSOL  Used  to  Solve  Constrained  Non-Linear  Programming  Problem 

□  Program  is  Fast.  Reliable,  Robust.  Flexible,  and  Readily  Extended 

□  Solution  Time:  Unperturbed  (sec.).  Perturbed  (min.) 


UcOonnmU  Dot$gki»  Spae»  Sy9t»m0  Company  •  Houaton  Dhrlahn  Paga22of22 
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MUUMP 


Multi-Impulse  Trajectory  and  Mass  Optimization  Program 


Daria  German 


Science  Applications  International  Corporation 


MULIMP  GENERAL  DESCRIPTION 


DESIGNED  TO  COMPUTE  A  MULTI-TARGETED  TRAJECTORY  AS 
A  SEQUENCE  OF  "TWO-BODY"  SUBARCS  IN  A  CENTRAL 
GRAVITATIONAL  FIELD  USING  KEPLER  AND  LAMBERT 
ANALYTICAL  SOLUTION  ALGORITHMS 


BODIES  MAY  BE  PLANETS  (ORBITAL  ELEMENTS  STORED 
INTERNALLY).  ASTEROIDS  OR  COMETS  (ORBITAL  ELEMENTS 
CONTAINED  IN  ASTCOM.ELM  RLE).  OR  FICTITIOUS  (ELEMENTS 
INPUT  BY  USER) 


CENTRAL  BODY  MAY  BE  THE  SUN.  ANY  OF  THE  9  PLANETS  OR 
AN  ARBITRARY  BODY  (GRAVITATIONAL  CONSTANT  INPUT  BY 
USER) 


UP  TO  19  SUBARCS  MAY  BE  SPECIFIED 


VARIABLES  IN  OPTIMIZATION  SEARCH 


•  TIMES  (DATES)  OF  THE  NODAL  POINTS 
CONNECTING  TRAJECTORY  SUBARCS 


•  POSITION  COORDINATES  OF  MIDCOURSE  AV 
POINTS  NOT  MADE  AT  AN  EPHEMERIS  BODY 


DEPARTURE  CONDITIONS 


•  RENDEZVOUS  DEPARTURE  IN  WHICH  CASE  THE  FIRST 
IMPULSE  AVi,  IS  EQUAL  TO  THE  HYPERBOLIC  EXCESS 
SPEED  Vhl 


•  PARKING  ORBIT  DEPARTURE  IN  WHICH  CASE  THE  FIRST 
IMPULSE  IS  THAT  NECESSARY  TO  ATTAIN  THE 
HYPERBOLIC  EXCESS  SPEED  FROM  THE  PARKING  ORBIT 
(fp.  e)  WITH  THE  MANEUVER  ASSUMED  TO  BE  COPLANAR 


•  A  -FREE*  DEPARTURE  IN  WHICH  CASE  THE  RRST  AV 
IMPULSE  IS  EXCLUDED  FROM  THE  PERFORMANCE  INDEX 


•  A  GRAVITY-ASSIST  DEPARTURE  IN  WHICH  CASE  THE 
APPROACH  HYPERBOLIC  VELOCITY  VECTOR  MUST  BE 
SPECIRED  BY  INPUT 


INTERMEDIATE  TARGET  CONDITIONS 

ARRIVAL 

•  RENDEZVOUS 

•  ORBIT  CAPTURE  (ORBIT  IS  USER  DEFINED) 

.  UNCONSTRAINED  FLYBY  SPEED 

.  CONSTRAINED  FLYBY  SPEED  (HYPERBOLIC  FLYBY) 
SPEED  IS  USER  INPUT 

DEPABTUBE 

•  RENDEZVOUS  DEPARTURE 

•  ORBfT  DEPARTURE  (ORBIT  IS  USER-DERNED) 

OTHER 

•  GRAVITY-ASSISTED  SWINQBY 


r 


GRAVITY-ASSISTED  SWINGBY 


•  MODEL  IS  FORMULATED  WITH  POWERED  MANEUVER  AS 
THE  GENERAL  CASE 


•  AV  WILL  OFTEN  ITERATE  TO  ZERO  VALUE  IF  THE  PROBLEM 
IS  NOT  OVERLY  CONSTRAINED  BY  SWINGBY  DATE  AND 
DISTANCE 


•  USER  OPTION  TO  SPECIFY  POWERED  MANEUVER  LOCATION 

•  INBOUND  ASYMPTOTE 

•  PERIAPSIS 

•  OUTBOUND  ASYMPTOTE 

•  BEST  CHOICE 


TERMINAL  TARGET  CONDITIONS 


•  RENDEZVOUS 

•  TARGET-BODY  ORBIT  CAPTURE 

•  SATELLITE  ORBIT  CAPTURE 

•  UNCONSTRAINED  FLYBY 

•  CONSTRAINED  FLYBY 

•  SPEaHED  ORBIT  ELEMENTS  (a.e.1)  RELATIVE  TO  CENTRAL 
BODY;  RNAL  TARGET  MUST  PROVIDE  GRAVITY-ASSIST 


16 


"FREE'  MIDCOURSE  AV  POINTS 


MIDCOURSE  VELOCITY  CHANGES  MAY  BE  MADE  AT 
INTERIOR  IMPULSE  POINTS  NOT  OCCURRING  AT  AN 
EPHEMERIS  BODY.  THESE  MIDCOURSE  AV  POINTS  MAY 
BE  INCLUDED  IN  TWO  WAYS: 


.  TIME  AND  POSITION  COORDINATES  MAY  BE 
ESTIMATED  AND  INPUT;  ON  USER  OPTION, 
THE  TIME  AND/OR  POSITION  COORDINATES 
WILL  BE  OPTIMIZED 


•  AUTOMATIC  IMPULSE  ADDITION  MAY  BE 
REQUESTED 


MULTIPLE  REVOLUTIONS 
AND 

RETROGRADE  MOTION 


MULTIPLE  REVOLUTIONS  AND/OR  RETROGRADE  MOTION  ARE 
SPECIFIED  BY  INPUT 


TWO  OPTIONS  ARE  PROVIDED  FOR  HANDLING  MULTIPLE 
REVOLUTION  SOLUTIONS: 


•  THE  NUMBER  OF  COMPLETE  REVOLUTIONS  AND 
ENERGY  CLASS  MAY  BE  SPECIFIED 


•  THE  MAXIMUM  NUMBER  OF  COMPLETE 
REVOLUTIONS  TO  CONSIDER  MAY  BE  ENTERED  IN 
WHICH  CASE  ALL  INCLUSIVE  SOLUTIONS  WILL  BE 
EXAMINED  AND  THE  'BEST  ONE  SELECTED  ON 
THE  BASIS  OF  A  VELOCITY  CHANGE  CRITERIA 


m&D  ENHANCEMENTS 


•  ADDITION  OF  THE  UNPOMEREO  SPECIFICATION  FOR  PLANETARY  SH1H6BYS 

•  A  NEW  DEPARTURE  OPTION  OF  SPACE  STATION  LAUNCHES 

•  A  NEW  TERMINAL  ARRIVAL  OPTION  OF  3-IMPULSE  PLANET  ORBIT  CAPTURE  TO  A 
FINAL  ORBIT  DETERMINED  BY  USER  INPUT  r^.  r..  AND  INCLINATION 

•  INCLUSION  OF  JPL  SATELLITE  EPHEMERIOES  ROUTINES  FOR  MOST  NATURAL 
SATELLITES 

•  CONVERSION  OF  THE  WORKING  COORDINATE  SYSTEM  FROM  EM050  TO  J2000 

•  ABILITY  TO  CONSTRAIN  TOTAL  TRANSIT  TINE 


RKSA  -  Applications  Branch 

Phillips  Laboratoiy 
Applications  of  POST 


James  B.  Eckmann 
SPARTA,  Inc. 
Phillips  Laboratory  SETA 
Edwards  AFB,  CA 

10  Aug  92 


RKSA  -  Applications  Branch 


Presentation  Overview 

•  Organization  and  Mission 

•  Simulation  Work  Environment 

•  Summary  of  POST  Models 

•  Applications  and  Some  Results 

•  Future  Plans 


RKSA  -  Applications  Branch 


Organizational  Hierarchy 


•jvmtAmotmcs 


AM.  CA 


RKSA  -  Applications  Branch 


Propulsion  Directorate  (RK)  Mission 

•  Provide  propulsion  technology  and  expertise  for  U.S.  space  and  missile  systems. 

•  Be  a  center  of  excellence  in  propulsion  research  and  development. 

•  Develop  a  broad,  advanced  technology  base  for  future  propulsion  system  designers. 

•  Demonstrate  propulsion  concepts  for  current  systems  designers. 

•  Assist  in  solving  operational  problems. 


9^ 

RKSA  -  Applications  Branch 

System  Support  Division, 
Applications  Branch  (RKSA) 

•  Mr.  Raymond  Moszee,  Branch  Chief 

•  Capt.  Tim  Middendorf 

•  1  Lt.  Paul  Castro 

•  2  Lt.  Naftali  Dratman 

•  Mr.  Francis  McDougall 

•  Mr.  Gerry  Sayles 

•  Ms.  Pamela  Tanck,  SPARTA 

•  Mr.  James  Eckmann,  SPARTA 

•  Maj.  Leo  Matuszak,  AF  Reservest 


RKSA  -  Applications  Branch 


RKSA  Simulation  Environment 


•  Integrated  Tools 

•  Ethernet  Network 

(TCP/IP) 

•  Connectivity  Software 

(NFS,  Versaterm  Pro) 


RKSA  -  Applications  Branch 


Integrated  Analysis 


Persuasion 


♦  POST  ♦ 


Sizing,  Geometry 
Structures  and  Weights 
Presentation 


MODEC 


TTTTrm 


APAS 


CONSE 


EVA 


SMART 
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RKSA  -  Applications  Branch 


Summary  of  POST  Models 

•  Atlas  II 

•  Delta 

•  Titan  IV  * 

•  Space  Shuttle  * 

•  Pegasus  * 

•  Ground  Based  Interceptors 

•  Small  ICBM 

•  Minuteman  III 

•  National  Launch  System  (NLS) 

•  Single  Stage  To  Orbit  (SSTO) 

•  Delta  Clipper  -  VTVL  Concept  ♦ 
•RASV- HTML  Concept 

•  RST  -  VTHL  Concept  * 


RKSA  -  Applications  Branch 


Applications  of  POST  Models 


•Alias  II 

•  Monopropellants:  High  Energy  Density  Mater  (HEDM)  propellants;  Composite  shroud 
•DelU 

•Titan  IV 

•  Soviet  RD-I70  straps  LRB's  to  replace  SRB's 

•  Space  Shuttle 

•  Qean  propellant  SRM's 

•  Pegasus 

•  Advanced  Liquid  Axial  Suge  (ALAS)  u  4th  suge;  Potential  booster  for  NASP  program 
flight  test  experiment 

•  Ground  Based  Interceptors 

•  Single-stage,  two-suge,  three-suge,  and  dual-pulse  motor  boosters;  Standard  Missile  and 
SRAM  2  biMsien  for  l^AP  testt 

•SICBM 

•  Advanced  ICBM  studies  baseline 
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RKSA  -  Applications  Branch 


Minuteman  HI  Model 


APPLICATION: 


*  The  ICBM  system  of  the  future.  Baseline  for  assessing  advanced 
technology  payoffs. 

•  Clean  Propellant  Trade  Studies 

•  Impact  of  Reducing  the  Number  of  Warheads 

•  Two-stage  Missile  Studies 
CONSTRUCTION; 

*  Objective  Function:  Maximum  range  for  fixed  payload  or  Maximum 
payload  for  Hxed  range 

*  Constraints  (2):  Maximum  dynamic  pressure.  Minimum  re-entry  angle 

*  Control  Variables  (6):  Pitch  rate  at  motor  ignition  for  each  of  3  stages; 
Tune  at  which  inertial  attitude  is  held  constant  for  each  of  the  3  stages 

*  Phasing  Events  (16):  3  motor  firings,  3  stage  seperations,  initial  pitch 
over,  3  constant  attitude  segments,  3  ballistic  flight  segments,  payload 
shroud  seperation,  atmospheric  re-entry,  ground  impact. 


RKSA  -  Applications  Branch 


Sample  Minuteman  III  Results 


RKSA  -  Applications  Branch 


National  Launch  System  (NLS)  Vehicles 


Trade  Studies  Performed: 
Castor  120's  on  NLS-3 
Mixture  Ratio 
Tank  Sizing 
Thrust  Level 
Engine  Out 
Throttling  Effects 
Staging  Algorithms 
Upper  Stages 


A 


y  ay  ay  0001,1 


A 


NLSS 


NLtt 


tatt 


QLPW.  MM  fij 

iSKtt  ISIM  1tX» 


4jI 

ITVtl 


RKSA  -  Applications  Branch 


Single-Stage-To-Orbit  (SSTO)  Models 


McDonnell  Douglas  Vertical 
TakeofT  Vertical  Landing 
(VTVL)  Delta  Clipper  model 
developed  and  provided  by 
NASA/Langley 


Rockwell  Intemational  Vertical 
Takeoff  Horizontal  Landing 
(VTHL)  Reusable  Space  Transport 
(RST)  model  developed  by 
Rockwell  and  provided  by 
NASA/Langley 


Boeing  Horizontal  Takeoff 
Horizontal  Landing  (HTHL) 
Reusable  Aerodynamic  Space 
Vehicle  (RASV)  model 
develop^]  in-house 


RKSA  -  Applications  Branch 


Future  Plans 

•  DEVELOPE  A  COMPLETE  VEHICLE  SIMULATION  CAPABILITY 

•  Apply  SMART  and  CONSIZE  to  current  analysis  tasks 

•  Complete  integration  of  Silicon  Graphics  machines 

•  Develop  a  cost  analysis  capability 

•  Continually  evaluate  new  analysis  tools 
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OTIS  Advances  at  the  Boeing  Company 
Steve  Paris 

Boeing  Defense  &  Space  Group 


Boeing 
Dotonso  S 

Space  Group 


Optimal  Trajectories  by  Implicit  Simulation 
(OTIS)  Development 


Boeing 
Defense  B 
Space  Group 


OTIS  Modes 


Mode  4  Trajectory  Optimization 

Direct  Trajectory  Optimization 

Optimai  Controi  1 

using  Nonlinear  Programming 

and  Collocation 

^ - h=hf 

_ -1- 

Mode  3  Trajectory  Optimization 
Parameter 


Mode  2  Targeting 
Shooting 


Boeing 
Defense  i 
Space  Group 


Next  Wave  of  OTIS  Advances 


Goal 


4  fold  overall  reduction 

OTIS  runtimes  reduced  by  10  fold 


Boeing 
Defense  S 
Space  Group 


OTIS  3.0  Lunar  Test  Case 


Explicit  Trajectory  Generation  Optimize 

•  Launch  Date 


•  A  parking  orbit 
•t0(burn1) 


\ 


24  hour  Orbit 


Low  Earth  orbit  to  high  polar  Earth  Orbit  (24hr). 
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Boeing 
Defense  A 
Spece  Group 


OTIS  Elements 


Boeing 
Defense  A 
Space  Group 


OTIS  3.0  Provides  Extreme  Flexibility 


•  Global  Constraints 

•  Analytical  Arcs 

•  Phase  Dependent 

•  Equations  of  Motion  (EQM) 
-  Control  Variables 

•  Quadrature  Variables 


BOOST 

FNgMPMheOM 
meh  and  Yaw  Control 


COAST 


Analytie 
NQ  vonvois 


nigMPaltiEQM 

An^  of  Attack  A  Bank  Controls 
ItoSMoVselor 
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Boeing 
Defense  & 
Space  Group 


Future  Trends 


Problem  Set-Up 


Data  Conditioning 


Program  Execution 


Results  interpretation 


(X  -  Window  Interfac^ 


Off-the-Shelf  Software 

•  Transform  by  Spyglass 

•  Excel  by  Microsoft 


New  Code 

•  Object  Oriented 

•  Software  Packages 
Technoiogies 

•  Sparse  Matrix  Methods 

•  Defect  Formulation 

•  Singuiar  Arcs 

•  Node  Piacement 


Off-the-Sheif  Software 

•  AGPS  -  Ribbon  Plots 

•  Agiie_Vu  >  Animated  Trajectories 


Boeing 
Defense  & 
Space  Group 


Summary 


•  Boeing  Continues  OTIS  Development 

•  Focus  on  Speed  &  Usability 

•  Exploit  Off-the-Shelf  Software 

•  Goal  is  a  "Better"  OTIS 
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—  4FS'n  ■■■  11.1  11.-  —  .  ■■■  .ini  - .  — M 

OTIS  ACTIVITIES 
AT 

MCDONNELL  DOUGLAS  SPACE  SYSTEMS  COMPANY 

R.  L  NELSON 
10  AUGUST  1992 

McDonnell  Douglas  Space  Systems  Company 

n.  N«Mon/3 
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OTIS  ACTIVITIES 


•  Applications 

•  OTIS  Project  Development  (PD  1-301)  at  MDSSC-HB 

•  Launch  Vehicle  Sizing:  ELVI8/OTIS 

•  OTIS  upgrades  for  Wright  Labs-AFB 


McDonnell  Douglas  SpacB  Systems  Company 

A  N«lMn/2 


ADVANCED  APPLICATIONS 
PD  1-301  OPTIMIZATION  TECHNIQUES  FOR 
ADVANCED  SPACE  MISSIONS 

■  MDSSC-HB  II  — — 

□  Theater  High  Altitude  Area  Defense  (IXAAD) 

□  ENDO/EXO  Atmospheric  Interceptor  (E^i) 

□  HEDI 

□  DELTA 

Q  National  Aerospace  Plane  (NASP) 

Q  SSRT 
Q  Asrobrakss 

□  Hypersonic  Advanced  Weapon  (HAW) 

□  Fighter  Aircraft 

e  Evasive  maneuvers 
e  Agility 

□  Military  Space 

Q  Space  Transfer  Vehicles 
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P0 1-301 
St3mOCKV2m.N* 


TASK  FLOW 

PD  1-301  OPTIMIZATION  TECHNIQUES  FOR 
ADVANCED  SPACE  MISSIONS 

MDSSC^HB 


TASK  1 -APPROACH  (1 992-1 993)-NLP  THEORY 
PD  1-301  OPTIMIZATION  TECHNIQUES  FOR 
ADVANCED  SPACE  MISSIONS 

M  MDSSC-HB 

□  D«v«lop  Strong  robust  globslly  convsrgent  nonlinear  optimizers 

•  Dense  and  sparse  optimizers 

•  NZ90L,  a  dense  optimizer 

•  Initial  feasibility  algorithms 

•  Min  •  Max  optimizer 

e  NZSPARSE,  a  sparse  optimizer 

e  Dr.  Philip  E.  Gill,  Professor,  University  of  California,  San  Diego 
Dr.  Michael  Saunders,  Research  Professor,  Stanford  University 


PD 1-301 
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TASK  1  -  PROGRESS-NLP  THEORY 
PD  1-301  OPTIMIZATION  TECHNIQUES  FOR 
ADVANCED  SPACE  MISSIONS 

MDSSC-HB 

□  D*v»lopMl  State  of  the  Art  optimizer,  NZSOL  ( dense ) 

•  NPSOL2.1 

•  NPaOL4.02 

•  Continual  teeting 

e  Tuned  NZSOL  for  OTIS  type  problems 

□  BREAKTHROUGH  ALGORITHM:  NZSPARSE  ( sparse ) 

•  Theoretical  formulation 

•  Development  and  checkout 
e  MINOS 

Q  Modified  OTIS  structure  to  accept  dense  /  sparse  optimizers 


P0 1-301 

TASK  2  -  APPROACH  (1 992-1 993)-NUMERICAL  ALGORITHM 
PD  1-301  OPTIMIZATION  TECHNIQUES  FOR 
ADVANCED  SPACE  MISSIONS 

MDSSC-HB 

Q  Algorithms  for  OTIS 
e  Automatic  scaling 

e  Automatic  node  placement  ( University  of  Illinois) 
e  Automatic  tabular  data  smoothing 
e  Lagrange  multiplier  Interpretation  (Continuous  /  discrete) 
e  Minimum  curvature  cubic  control  splines 
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TASK  2  -  PROGRESS  -  ALGORITHMS 

^AFS  —  !■  ■■■■■  ■'  —  -■ 


•  Tabular  Data  Smoothing 

•  Enhancad  Velocity  Lots  Model  for  Launch  Vehicles 
e  Generalized  Stage  -  Phase  Concept  for  Sizing 

e  Automatic  Node  Placement 


TASK  3  -  APPROACH  (1 992-1 993)-LOW  THRUST 
PD  1-301  OPTIMIZATION  TECHNIQUES  FOR 
ADVANCED  SPACE  MISSIONS 

.  MD$SC-HB  HIM . .  ■■■  !■!  . . . . . .  ■■■!  ■  111 

□  Develop  restricted  and  gsneral  3-body  equations  of  motion 

□  Boundary  conditions  and  coordinate  systems 

□  Quantify  the  transition  region  for  earth-moon  low  thrust  /  weight 
transfers 

□  SECKSPOT/ NASA  Code -COSMIC  Library 
e  Strong  gravity  field 

e  Orbit  averaging  techniques 
Q  QT2  interplanetary  code 
Q  Dr.  Richard  Shi,  MDSSC-HB 


PD 1-301 
asamocKvimjw- 
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TASK  3  -  PROGRESS-LOW  THRUST 
PD  1-301  OPTIMIZATION  TECHNIQUES  FOR 
ADVANCED  SPACE  MISSIONS 

MDSSC-HB  ■■ 


□  K«y  IdM  to  Mlv*  probitm 

•  Tho  oxlotonco  of  tho  JooobI  (onorgy)  integral  for  the  restricted 
three-body  problem  will  eld  us  In  the  general  three-body 
problem. 

•  Zero  velocity  or  zero  energy  curves  are  the  regions  where 
the  low  thrust  earth-moon  transfers  are  possible 

e  No  integral  available  for  the  general  3-body  problem 

Q  Develop  for  OTIS 

e  3-body  equations  of  motion 

e  boundary  conditions 

e  coordinate  systems 


P0 1-301 

SSWROCKV2mLH«n>’ 


PD  1-301  OPTIMIZATION  TECHNIQUES  FOR 
ADVANCED  SPACE  MISSIONS 


VKB1S63 


Moon 


Earth-Moon  Transfer 
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AFS 


ELVIS  /  OTIS 
ARCHITECTURE 


Atm— wq 


UeDomtttt  Douglat  Spsc»  Sytlama  Comptny 
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OTIS  UPGRADES  FOR  WRIGHT  LABS 


—  AFS 


•  TASK1:  NZSOL 

•  TASK  2:  Automatic  Variable  Scaling 

a  TASKS:  Variable  Names  for  NZSOL  Output 
a  TASK  4:  Minny  Heating  Model 


McDoniwU  Douglas  Space  Syatems  Company  ^ 

a  N*iMfl/a 
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Advances  in  Trajectory 

Optimization  Using  Coiiocation 
and  Noniinear  Programming 


Bruce  A.  Conway 

Dept,  of  Aeronautical  &  Astronautical 

Engineering 

University  of  Illinois 

Urbana,  IL 


August  1992 
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Outline 


Introduction 

Progress  to  Date  •  Theory 
Progress  to  Date  *  Solved  Problems 
Continuing  and  Proposed  Research  -  Problems 


Progress  -  Theory 


1 .  Use  of  costates  to  improve  an  qstim^  trajectory. 

Lagrange  multipliers  for  the  discrete  (NPSOL)  solution  are  a 
representation  of  the  Lagrange  multipliers  of  the  continuous  case, 
(bright  &  Conway,  JGC&D  15,  No.  4, 1992) 

Knowledge  of  the  Lagrange  multipliers  allows  a  posteriori  determination 
of  the  optimality  of  the  solution,  e.g.,  can  examine  the  switching  function. 

2.  Generalized  defects 

Can  be  used  when  the  differential  equation  for  a  state  variable  is 
integrable,  e.g.,  on  a  coast  arc. 

May  significantly  reduce  the  number  of  NLP  parameters  and  hence 
execution  time. 

3.  Coordinate  transformation  within  the  H-P  structure 

Necessary  for  orbit  transfer  when  changing  sphere  of  influence 
Keeps  state  variables  near  one  order  of  magnitude,  as  NPSOL  prefers 
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4.  Method  of  parallel  shooting 

Replaces  single  Hermite-Simpson  "integration  step”  with  multiple 
Runge-Kutta  steps  allowing  use  of  larger  intervals. 

Results  in  smaller  NLP  problems  for  a  given  accuracy. 

5.  Automatic  node  placement 

Computer  solves  a  succession  of  NLP  problems  in  which  additional 
nodes  are  inserted  as  needed  to  acheive  a  given  accuracy. 

More  efficient  than  using  a  uniform  distribution  of  nodes 

6.  Neighboring  optimal  feedback  control 

Determines  gains  for  linear  feedback  controller  to  yield 
neighboring  optimal  controller 

Unnecessary  to  solve  NLP  problem  for  small  change  in  initial 
or  terminal  conditions 

Feedback  gain  history  easily  loaded  into  small  memory 


Illustration  of  Generalized  Defects 
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Illustration  of  Parallel  Shooting 


Progress  -  Solved  Problems 

1.  optimal  iow-thrust  escape  trsyectory  (Enright  Ph.  D.  thesis) 

2.  Optimal  2  and  3  burn  circle-circie  low-thrust  rendezvous 
(Enright  Ph.  0.  thesis) 

3.  Optimal  low-thrust  Earth-Moon  transfer  (Enright  Ph.  0.  thesis) 

4.  Optimal  spacecraft  detumbling  (A.  Herman  M.S.  thesis) 

5.  Optimal  low-thrust  insertion-  Mars  Observer  (Enright  Ph.  0.  thesis) 

6.  Optimal  20  and  30  direct  ascent  time-bounded  interception 
(J.  Oowney  Ph.  0.  thesis) 

7.  Neighboring  optimal  feedback  control  for  continuous-thrust  ascent 
maximizing  horizontal  velocity  (F.  Chen  Ph.  0.  research) 
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Low-1  hrust  Minimum  Fuei  Escape 


Mflthod  _ Vari^ales - CEU 


Harmila/SunpMn  (60) 

427 

190  sac 

ParalM  ahooiing  (34  x  3) 

385 

96  sac 

ParaHal  shooting  (5  x  20) 

270 

72aae 

ri=  1 

ai«  0.0125 
=  16ji 

All  in  canonical  units 


Optimai  2  and  3  Burn  Circle-Circie  Rendezvous 
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COPY  AVAILABLE  TO  DTIC  DOES  NOT  PERMIT  FULLT  LEGIBLE  REPROi/LwilON 


thnist  auigle  (degrees) 


Optimal  Low-Thrust  Earth-Moon  Transfer 


rig.  f  Cant  ■■■■  iTMifM 


Optimal  Low-Thrust  Earth-Moon  Transfer,  cont'd 
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optimal  Spacecraft  Passivation  (Detumbling) 


View  of  OMV  /  Disabled  Satellite  System 


Spacecraft  Passivation,  cont'd. 


Results  from  the  * 

I 

TPBVP  solver  1 

» 

1 


Results  from  the 
NLP  method 


External  Torque  Histories 

A9 

COPY  AVAILABLE  TO  DTIC  DOES  NOT  PERMIT  FULLT  LEGIBLE  REPRODLCTION 


Optimal  2D  &  3D  Direct  -Ascent  Interception 


•  Target  is  assumed  to  be  in  a  general  Keplerian  orbit  with 
orbital  elements 


£■*^  =  [3,6,  i,  O,  to,  f] 


•  Geometry  of  the  problem 


Continuing  Research  -  Problems 


1.  Automatic  node  placement.  (A.  Herman) 

2.  Optimal  very-low-thrust  trajectories  (W.  Scheel) 

3.  Optimal  Earth-Mars  low-thrust  transfer  including  escape  and 
arrival  spirals  and  coordinate  transformations  at  sphere  of 
influence  of  each  planet.  (S.  Tang) 

4.  Neighboring  optimal  feedback  control  for  complex  problems 
Automation  of  NOFC  using  symbolic  programming  (F.  Chen) 

5.  Optimal  trajectories  for  interception  of  Earth-crossing  asteroids 
(B.  Conway) 
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Accuracy  of  OTIS  Solutions 


OTIS  :  Optimal  control  solutions  via  direct  transcription 

Combination  of  collocation  and  nonlinear  programming 
Question : 

How  do  OTIS  solutions  compare  to  the  "exact”  or  TPBVP 
solutions  ? 

•  Analytical  Approach 

estimate  adjoint  variables,  examine  discretized  necessary 
conditions 

•  Engineering  Approach 

numerical  comparison  of  collocation  solution  to  exact 
solution  for  a  representative  problem 
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FLIGHT  PATH  OPTIMIZATION  OF 

AEROSPACE  VEHICLES  USING 

OTIS 

Rajiv  S.  Chowdhry 

Lockheed  Engineering  &  Sciences  Company 

MS  489 

Aircraft  Guidance  &  Control  Branch 

NASA,  LaRC,  Hampton  VA. 
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Outline 


•  Accuracy  of  OTIS  solutions. 

•  Overview  of  OTIS  applications  at  AGCB 


Accuracy  of  OTIS  Solutions 

OTIS  :  Optimal  control  solutions  via  direct  transcription 

Combination  of  collocation  and  nonlinear  programming 
Question : 

How  do  OTIS  solutions  compare  to  the  "exact”  or  TPBVP 
solutions  ? 

•  Analytical  Approach 

estimate  adjoint  variables,  examine  discretized  necessary 
conditions 

•  Engineering  Approach 

numerical  comparison  of  collocation  solution  to  exact 
solution  for  a  representative  problem 
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Example  :  ALS  Ascent  to  Orbit 


Example  Problem : 

Steer  a  two  stage  launch  vehicle  from  a  given  Initial  condition  to  a 
specified  target  orbit  In  minimum  fuel. 


•  Exact  or  TPBVP  solution  available  in  literature  (  Ref.  Hans 
Seyv^d  and  E.  M.  Cliff ) 

•  Care  was  taken  to  keep  the  vehicle/atmosphere/planet  models 
same  in  OTIS 

•  Only  solution  methodologies  were  different 


t|  MC 


»(tf)Kgs 


V«lodty(l,) 

m/Mc 


AIWuda  ( If )  km 


CPU  Tim«  (i 
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32  nodM 

40  nodM 

477.2 
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149.991 

7995 

7955 
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149.99 

149.74 

147.90 

147.99 

274.59 

274.11 
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149.53 

149.99 

147.9 

147.9 

74 

377 

935 

1975 
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477.2 


149.900 


7957 


149.2 


277.91 
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Figure  [  1  ].  Comparison  of  optimal  ALS  ascent  with  OTIS 
solution,  mass  (kg)  vs.  time. 


Figure  {2].  Comparison  of  optimal  ALS  ascent  with  OTIS 
solution,  altitude  above  spherical  Earth  vs.  time. 


5^ 


Figure  [3].  Comparison  of  optimal  ALS  ascent  with  OTIS 
solution.  Earth  relative  velocity  (m/sec)  vs.  time. 


Figure  [4].  Comparison  of  optimal  ALS  ascent  with  OTIS 
solution,  local  horizontal  flight  path  angle  (deg)  vs.  time. 
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•  a  T - 1 - 1 - 1  I  r 

0  100  200  300  400  500 

Time  (sac) 


Figure  [S].  Comparison  of  optimal  ALS  ascent  with  OTIS 
solution,  thrust  vector  angle  (deg)  vs.  time. 


OTIS  Applications  at  AGCB 

•  Fuel  efficient  ascent  for  SSTO  airbreathing  hypersonic  vehicle. 

•  fuel  optimal  path  definition  for  G&C  studies 

•  HL-20  abort  maneuvers :  ELSA  ( Efficient  Launch  Site  Abort ) 

•  Parameter  sensitivity  studies  to  support  design  activities. 

•  Guidance  algorithm  development  &  real  time  validation. 

•  ALS  ascent  for  OHS  calibration. 

•  Optimal  maneuvers  for  a  high  performance  fighter  aircraft  (HARV) 
in  air  combat  situation. 
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Conclusions 


For  the  ALS  Ascent  Problem  : 

•  Excellent  match  of  the  collocation  solution  to  the  TPBVP 
solution 

•  Relatively  quick  turnaround  time  for  OTIS  solutions 

•  Very  robust  to  initial  guesses 
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Trajectory  Optimization  of  Launch  Vehicies 
at  LeRC:  Present  and  Future 


Presented  by 
Koorosh  Mirfakhraie 


at 

Workshop  on  Trajectory  Optimization  Methods  and  Applications 

Hilton  Head.  SC 
August  10,  1992 


AN  ALEX 

CORPORATION 


KM  i/iwn 


1V/\SA  Lewis  Research  Center 
Advanced  Space  Analysis  Office 
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Outline 


Introduction 

Present  method  of  solution  and  code 
Capabilities  of  the  present  code 
Motivation  for  replacing  the  code  (and  method) 
Examination  of  methods  using  collocation 


Introduction 


Trajectory  optimization*  of  ELV's  at  the  Advanced  Space 
Analysis  Office  at  LeRC  is  performed  for: 

•  Mission  design  for  approved  programs 

•  Feasibility  and  planning  studies 

•  Corroboration  of  contractors’  data  for  NASA  missions 
flown  on  Atlas  and  Titan 


*  Trajectory  optimization:  Maximizing  the  final  payload 
subject  to  a  set  of  intermediate  and  final  constraints. 


ANALEX  IMASA  LawJ*  Research  Center 

awoRATKJN  Advanced  Space  Analysis  Office 

KM  VWfl 
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Introduction  (Cont’d) 


Mission  profiles  for  launch  vehicle  systems  with  booster  and 
upper  stages  include: 

•  Launches  from  ER  and  WR 

•  LEO,  GTO,  and  GSO  insertion 

•  Interplanetary  escape  trajectories 

•  Orbit  transfers 


Present  Method  of  Solution  and  Code 


Calculus  of  Variations  is  used  to  formulate  the  problem. 
The  resulting  two  point  boundary  value  problem  is  solved 
using  a  Newton-Raphson  algorithm. 

The  computer  program  (DUKSUP)  was  written  entirely  at 
LeRC  during  1960’s  and  early  70’s. 

DUKSUP  is  a  3-D.O.F.  code  written  for  performance 
analysis  of  multi-stage  high-thrust  launch  vehicles. 


ANALEX 

CORPORATION 

KM  t/IVtl 


Lewis  Research  Center 
Advanced  Space  Analysis  Office 
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DUKSUP  Features 


Detailed  modeling  (e.g.,  propulsion  and  aerodynamic) 
of  a  launch  vehicle  is  possible. 

A  variety  of  constraints  can  be  imposed  on  the  model. 
They  include: 

-  Instantaneous  and  total  aerodynamic  heating 

-  Maximum  dynamic  pressure 

-  Parking  orbit  parameters  (e.g.,  radius  of  perigee, 
energy,  velocity,  etc.) 

-  G-iimit  staging 


DUKSUP  Features  (Cont’d) 


Several  in-plane  and  out-of-plane  final  target  conditions 
can  be  specified  (e.g.,  energy,  radius,  true  anomaly, 
inclination,  declination  of  outgoing  asymptote,  etc.). 

Variables  free  for  optimization  include: 

-  Upper  stage  burn  and  coast  times 

-  ‘Kick  angle’ 

-  Payload  fairing  jettison  time 

-  Thrust  angle  in  the  non-atmospheric  flight 


ANALEX 


CORPORATION 


KM  vttm 


fUASA  Lewis  Research  Center 
Advanced  Space  Analysis  Office 
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Motivation  for  Replacing  DUKSUP 

Sensitivity  to  initial  guesses 

Difficulty  in  reformulating  the  C.OA.  problem  when 
adding  new  features  and  constraints  to  the  code 

Difficulty  in  modifying  and  expanding  the  code  due  to 
lack  of  documentaion  and  outdated  programming 
practices 


Examination  of  Methods  Using  Collocation 

•  Two  main  features  of  collocation  making  it  attractive  are 

•  Lack  of  sensitivity  to  initial  guesses 

•  Relative  ease  of  formulation 

•  Concerns  about  using  collocation  for  ELV  optimization  are 

-  Ability  to  handle  complex  modeling  requirements  and 
constraints  typical  of  ELV  flight 

-  Computer  run  time 

•  Fidelity  of  the  solution  vis  a  vis  C.O.V. 

•  Evaluation  of  collocation  uses  DUKSUP  as  the  benchmark 
for  comparison. 

IVI/NSA  Lewis  Research  Center 

KM  immi  Advanced  Space  Analysis  Office 
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Using  Collocation  (Cont’d) 


Available  collocation  codes  are  used  as  testbeds  with 
necessary  modifications. 

A  simple  LV  model  is  used  first  and  moved  progressively 
to  a  full  DUKSUP  model. 

Enright’s  orbit  transfer  program  was  used  for  the  first 
simple  model  comparison.  Results  matched  those  of 
DUKSUP. 

OTIS  is  used  for  the  more  sophisticated 
comparisons. 

OTIS  is  currently  used  to  model  an  Atlas  ll/Centaur  to 
LEO. 


ANALEX 


CORPORATION 


KM  viwn 


Lmvls  Research  Center 
Advanced  Space  Analysis  Office 
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Collocation  Methods  in  Regular  Perturbation  Analysis 
of  Optimal  Control  Problems* 


August  10«  1992 


Prepared  for 


Workshop  on  Trajectory  Optimization  Methods  and  Applications 
AIAA  Guidance,  Navigation,  and  Control  Conference 
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Anthony  J.  Calise**  &  Martin  S.K.  Leung 
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Overview 


Motivation 

Regular  Perturbation  Analysis 
The  Method  of  Collocation 

Hybrid  Collocation  /  Regular  Perturbation  Analysis  Approach 
Examples 

Duffing  Equation 

Launch  Vehicle  Guidance  Application  (presented  at  1*GNC>1) 


Motivation 


AMlyttetf  McUMtft  N— Brirri  M»lho4f 
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Advantages  /  Disadvantages 


Analytical  M^bods 

Approximates  solution  by  cxpanskm  in  an  asymptotic  series  in  a  snutt 
parameter 

Zero  order  problem  is  simpler  to  scdve  =>  Insight 
Higher  order  problems  are  linear 

Zero  order  problem  must  reasonably  approximate  the  full  order  problem 

For  practical  applications,  xero  order  problem  must  be  analytically 
tractable  or  reducible  to  a  simple  algebraic  problem 

Significant  amount  of  analysis  is  required  for  each  problem  formulation  of 
interest 


Advantages  /  Disadvantages  (continued) 

Numerical  (Collocation)  Methods 

Finite  element  method  that  enforces  interpolatory  constraints  at  specific 
points  within  each  element 

Simple  to  use  for  a  wide  variety  of  optimization  problems 
Large  dimensional  nonlinear  programming  problem 
No  general  guarantee  of  convergence 


Advantages  of  analytical  and  numerical  methods  are  fai  many  respects 
complimentary  in  the  sense  that  if  the  advantages  can  be  combing  in 
some  way,  then  most  of  the  important  disadvantages  for  real*tinie 
applications  can  be  removed. 
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Regular  Perturbations  in  OptinuU  Control 

Givco: 

dx/dt  8  rCx^t)  •»>  e  g(x,u4);  x(t«) »  x« 

j8(|>(X,t)lj^ 

Find  the  cootrol  that  minimizes  J  subject  to  the  terminal  time  constraints: 

V(x.t)lg^80 

Optimality  mndition; 

Hh  8  0  assuming  Hvu  >0  ^  u  s  U(x,  X,  t) 

where: 

H8XT{f+eg);  H(tf)  =  .(I>||^;  + 

tUydts.Hi;  Mtr)  =  Oxlt^ 


Regular  Perturbation  Analysis 

Based  on  a  simplified  modei  (when  e  is  set  to  zero) 

•  Treat  neglected  dynamics  as  perturbation 
-  Define  a  normalized  independent  variable*  Ts(t*t«)/T 

whereTstfto 

•  Compute  zero  order  solution 

Consider  an  asymptotic  series  in  x,  X*  and  T 

Evaluate  high  order  corrections  from  sets  of  nonhomogeneous* 
time-varying  linear  O.  D.  E's. 

dtL^kJ  L^ll  A22JL^|iJ  T*  .C2J  [Pik 

enforcing  all  boundary  conditions  to  lith  order 

Compute  feedback  control  at  current  time  (to)  using  x(to)  and  kth  order 
approximation  for  X(to) 
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Regular  Perturbation  Analysis  (continued) 


•  A's  and  C's  depend  only  on  the  zero  order  (k  s  0)  values. 

•  C's  are  the  explicit  correction  term  for  free  final  time,  T. 

•  Ps  are  the  forcing  functions  involvii^  lower  order  (k-1,  ^  1, 0)  terms. 

Solution: 


Higher  order  correction  involves  simple  operations  of  quadrature  and  soluthm 
of  linear  algebraic  equations 

Can  be  easily  modified  to  account  for  discontinuous  dynamics 


Solution  of  Optimal  Control  Problems  by  Collocation 


Methodology 

•  a  Unite  element  approach 

•  approximates  the  solution  with  interpolating  functions 
.  consider  first  order  polynomials 

x(t)  =  x^_i+Pj(t-tj_i)  ;  A,(t)  =  Xj_i +qj(t-tj_i)  ;J  =  1,  2,  ..,N 


•  enforce  the  derivative  constraints  at  the  mid  point  of  each  element 

i=(ij+ij_l)/2;  x=(»j+*j_i)/2;  2.*(Xj+X^|)/2 
i=(ij  +» j-i  /2;  *=(* j +» j_|  )/2;  X=(2,j  +X j_|  )/2 


•  N  is  the  number  of  elements,  xj  jod  Xj  are  nodal  values 

•  control  assumed  to  be  eliminatec!  using  optimality  condition 
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Hybrid  Collocation  /  Regular  Perturbation 

A  Regular  Perturbation  Formulation 

-  rewrite  the  actual  dynamics  as 

,dH  ,  .•  ,  dH  , 

x  =  Pj+e(^-Pj)  ;X*qj+e(~-qj) 

-  perturbathm  terms  are  zero  at  mid  point  of  each  element. 

•  for  cases  that  control  cannot  be  eliminated  explicitly,  use  an 
analytic  portion  n(x,  X,  u) 

o=n+e(~-n) 

du 

Carry  out  a  Regular  Perturbation  Analysis 

.  expand  about  the  zero  order  solution  (derived  from  collocation) 

•  provides  higher  order  corrections  to  collocation  solution 

•  Anther  exploitation  of  the  analytically  tractable  portion  of  the  dynamics 
will  result  in  more  intelligent  interpolating  functions  (see  simple  example) 


A  Simple  Example 


Dufllng’s  equation  in  first  order  form: 

i  =  v  ;x(0)  =  Xp 

v  =  -x-ax^+u  ;v(0)=Vp 

J  s  S*  xi(tf)  +  Sy  v2(tf)  +  +  uV2)  dt 

Notes: 

•  hardening  effect  is  given  by  the  nonlinear  term,  ax^ 

•  the  optimal  control  problem  ir.  a  fourth  order  example 

•  will  demonstrate  different  levtils  of  intelligent  interpolating  functions 
that  enhance  the  approximation  with  fewer  number  of  elements 
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Simple  Example  (continued) 


Level  0  Formulation; 

•  degeiicratecase,usesoiily  regular  perturbation  With  a  completely 
analytic  zero  order  solution 

•  let  e  s  a  s  0.4,  and  treat  the  nonlinear  terms  as  perturbations 

•  Si  s  Sv  s  100 


i  =  v 

♦  =  -x+u-ex^ 

Xy  =  — Xj 

Hy  =0  + Xy  -0 


X(0)=Xy 
V(0)=Vy 
Xy(t,)  =  2S,x(tr) 

X^(t,)  =  2S»v(tr) 

{H  =  X,v+Xy(-x  +  u-ex^)+l+u*/2}|tf  =0 


•  zero  order  problem  is  linear  and  time«invariant 

•  compute  up  to  second  order  corrected  solutions  (Fig's.  4.1  and  4.2) 

•  series  not  convergent,  most  accurate  approximation  is  first  order 

.  nonlinear  term  ax^  is  too  large  to  be  neglected  in  the  zero  order  problem 
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Simple  Example  (continued) 


UydiEarmuiiiiflfl; 

•  luie  hybrid  approach,  approximate  all  state  and  costates  as  piecewise 
linear  functions 


=  J  v,(t)=v,j_i  +  pyj(t-tj_,) ;  j  =  l,  ..,N 


•  number  of  unknowns  is  4N  5 

•  Istandlnd  order  corrections  are  computed  for  Nb3  (Fig's.  4i9- 4.12) 

•  discontinuity  in  slope  is  smoothed  as  order  of  correction  increases 

•  correction  by  regular  perturbatioo  analysis  allows  use  of  crude  number 
of  element  representation  in  the  zero  order  collocation  solution 
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Figure  4.9.  Level  1  Higher  Order  Results 
In  X  for  N  s  3. 


Figure  4.10.  Level  1  Higher  Order  Results 
in  V  for  N  K  3. 


Simple  Example  (continued) 


Level  2  Formulation: 

•  enhanced  level  1  formulation  by  Interpolating  only  those  variables  that 
have  nonlinear  coupling 

•  decompose  the  dynamics  as: 

dx/dts  V 

dv/dt  s  pvj  +  e{-x  -Lv  -ax3  -pvj)  j  =  lA  -»N 

dXi/dt  s  -I-  £{  A«(l  ■¥  3ax2)  •  q^j} 

dXv/dt  s  ’Xx 

-  number  of  unknowns  is  2N  -i-  5 

•  both  zero  and  first  order  results  for  N  =  2  are  superior  than  the 
N  s  3  results  for  the  Level  1  formulation  (Fig's.  4.13  •  4.16) 
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Simple  Example  (continued) 

Level  3  Formulation; 

>  enhanced  level  2  fonnulation  by  fully  utilizing  analytically  tractable 
portion  of  the  necessary  conditions 

•  decompose  the  dynamics  as: 

i  =  v 

♦  =  -z-X,+Pyj+e(-az^-Pvj)  ;j=l,  2, N 

Xji  =  Aiy +(|jj +e(3aXiyZ 
ly  =  -Xg 

•  similar  to  Level  0  except  for  additional  unknown  constants  pyjyQ:^ 

•  use  piecewise  constant  terms  to  approximate  the  nonlinear  parts 

•  both  zero  and  first  order  results  for  N  =  1  are  superior  than  the 
Level  0  case  (F^i's.  4.17  •  4.20) 

•  I  t  and  3  casga  demonstran!  the  use  of  intelligent  interpolating 
ftinctions 
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Figure  4.17.  Level  3  Higher  Order  Results 
in  X  for  N  s  1. 


Figure  4.18.  Level  3  Higher  Order  Results 
in  V  for  N  s  1. 


Fifiiit4.I9.  l.evel3Hi|herOR>erRetuIainX(forN>l. 


Alternative  Implementations 


Repeat  zero  order  soiution  and  perform  quadratures  at  each  control  update 
interval 

Or 

Compute  zero  order  solution  and  quadratures  off  line,  and  store  for  in-flight 
use 


Improves  reliability  and  computational  efficiency  with  some  loss  in  accuracy 
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Summary 


Beoefits  of  Hybrid  Approach: 

SigpiflcanUy  improves  a  collocation  solution 

First  and  higher  order  corrections  are  obtained  by  quadrature 

Intelligent  interpolation  functions  obtained  by  retaining  as  much 
of  the  analytically  tractable  portion  of  the  solution  as  possible 

Possible  to  implement  the  control  solution  so  that  the  zero  order 
solution  and  quadratures  are  performed  once  ofMine  and  stored 

Significantly  improve  a  regular  perturbation  solutieii 

Retain  more  of  the  nonlinearities  in  the  zero  order  problem  by 
using  finite  elements  and  collocation  to  construct  an  improved 
zero  order  solution 


Important  implications  in  real>time  guidance  applicutions 
Computational  efficiency  and  reliability 


AUTOMATIC  SOLUTIONS  FOR  TAKE-OFF 
FROM  AIRCRAFT  CARRIERS 


Lloyd  H  Johnson 
AIR-5301 2D 
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SMtioll  NAfCMfOO  NAtCMMO 


I_^LJNCM  Ff|_I<3MT  P*«XM 


NACCOffOO 


OlAtUUNCH 


444  SHIDU/rENOANT  Ut'NCH  METHOD. 

ThU  amliod  of  laiMKfciH  «*<«*•(<  «•  ••  iMf**  «  Mf* 

ofb«»  M«ct  tht  now  fnr  iMHidi  Miliod  biCMif 

<«m  pMignpfc  44.1).  WHk  dii*  wuchod.  tilt 
•iicvtft  it  co«pM  M  tlw  cttapvii «»«  ftftinf  ky  of 
%  «ri«t  topt  Wi^  Of  ytiMtm.  TM  WMIt  it  “V"  fhtptd 
•04  ft^oittt  («o  «o«  ficciflgt  o«  tht  tivcttfi.  ohtrttt  «M 
pto4twi  ottdt  onlf  OM  fo«  finitf  oo  clw  tircrafi. 

Tlw  boUbtch  4t*fet  ■••4  wtili  fari4h/pt<i4tot  Itwodi  it 
v«fo  f«^.  chtio.  Of  Mttil  linkt.  Tlw  itlttit  tkmmm  it 
tidwr  rittf  or  a  ttotioo  Wf.  Tlw  Iwllfctffc  otttwMy  ttcwii* 
OHOc  pokm  oo  tlw  tittftfi  it  wtH  ifi  wM  tlw  4»ch  tiM  m> 
ladwt  w  tlw  boMback  4trfc  dm.  Stcitoo  Vlll  4tiCTibtt 
tfpkd  bol4btck  io4  rtitttt  intwMitt 

That  iiwtlio4  of  Itoocbing  rtqwiftt  owootf  Imkiif  of 
dw  bn4lt  or  |Wft4t<it  to4  Uw  buMktfti  by  tbt  ctwfolt 
4tcb  cftv  tfttr  ifw  tiicrtft  Iwt  btta  wtw4  into  pmt- 
tmt  oo  tbt  cttfwk.  Wbto  ctw  tircftfi  rttcbti  ibt  to4  of 
cbt  ctttyoh  gowtf  rail  tii4  tbt  too  fofcv  4ttt|r«,  tbt 
bfi4H  or  ptodtM  4roft  frow  tbt  tifcrtft  tow  finiiift 
•04  it  bfMf bt  to  t  «of  on  tbt  flight  4tcb  bf  tbt  bci4lt 
artttttf  tytttio. 

4.1  LAUMCMM6  IQUirmCMT. 


44.4  CENERAt. 


4.Id>.4  CATAPULT*  A  dirvtf^riro,  tittw  ITF* 
ttttfiitl  it  tMo4  oo  oU  carritrt.  Sctaw,  ptbt4  tbt 
tbi^'t  botitfft  to  0  ttritt  of  Itrgt  tttow  roctimt*  it 
ftiootod  ta44«alf  to  tW  itoocbi^  eogiot  to  4rra«  rwo 
pttMOt.  TW  pitfoot  tft  diftctiy  ctootcltd  to  iW 
ctiofoit  too  fitciog  (broogb  tbt  tlocitd  cyliodtr  ooUt. 
Rtirocoto  it  otcfpUtbtd  hf  t  ttptrttt  byAw 
pofotif  mrtcciag  togioo. 

4.14.1  DECK  EQUIPMENT.  TWcwtpob  It  twoipM 
wiib  t  ibofrit  wbicb  loom  tloog  tbt  cacapoli  tncb 
dofiiii  ibt  iMMcbi^  optttiioo  tod  friotwitt  dw  cttt* 
poWm  fetct  freto  dw  cntpoli  tufiot  lo  dw  tircfMi 
tbfoogb  dw  hoiidi  btf.  brMIt  or  pttdttt  A  rwop.  tt> 
otfod  to  ibt  ciapob  urn  fining,  totbiii  dw  •irtnh  to 
rofl  om  dw  tow  fiitiog.  For  bridit  Itii^td  tirmfi.  tbt 
boldbtct  4wb  dm  pnn  idtt  woowrotw  tocbtrtgt  poittt 
for  tW  bol4bwb  wd  nttttt  tiwwbly.  Tbh  dtch  dm 
it  loctttd  aft  of.  tod  oo  dw  ctowHiot  of  ibt  rntpob 
ifocb.  TbtbridM  trrttwr  tiopa  tod  mtiot  dw  bridte  or 
ptodwn  dwr  it  it  tbtd  from  tlw  tirrrtfi.  Fintlly.  tbt  dtcli 
fdfft  cooirol  ptod  providtt  tbt  priowry  cootrol  pom 
Aw  uptiinog  tbt  caupoH. 
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LTV  Aerotpac*  and  Oefanaa  Company 
Aircraft  DIvlaion 


PROGRAM  OVERVIEW 


THE  CATAPULT  LAUNCH  SIMULATION 
CONSISTS  OF  FIVE  PHASES 

•  STATIC  BALANCE 

•  HOLDBACK 

•  CATAPULT  STROKE 

•  DECK  RUN 

•  FLYAWAY 


LTV  Aerospaca  and  Oefanaa  Company 
Aircraft  Division 


program  overview 


THE  CATAPULT  LAUNCH  SIMULATION  INCLUDES: 

•  CATAPULT  FORCES 
a  HOLDBACK  FORCES 

•  HIGH  FIDELITY  LANDING  GEAR  MODEL 

•  AERODYNAMIC  DATA  AS  A  FUNCTION  OF 

-  ANGLE  OF  ATTACK  OR  LIFT  COEFFICIENT 

-  NOZZLE  DEFLECTION 

-  THRUST  COEFFICIENT  OR  NOZZLE  PRESSURE  RATIO  (NPR) 

-  FLAP  DEFLECTION 

-  PITCH  TRIM  SURFACE  DEFLECTION 

•  GENERIC  FLIGHT  CONTROL  AND  STABILITY  AUGMENTATION 
SYSTEM 

•  LONGITUDINAL  THRUST  VECTORING 


RMCOM 
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LTV  Atrospac*  and  Dafansa  Company 
Aircraft  Division 


PROGRAM  OPTIONS 


PROGRAM  OPTIONS 

•  AUTOMATIC  WIND  OVER  DECK  SOLUTION 

•  SOLUTION  TERMINATION 

•  POWERSETTING 

•  FLAP  DEFLECTION 

•  FLIGHT  CONTROL  AND  STABILITY  AUGMENTATION  SYSTEM 

•  THRUST  VECTORING  CONTROL  SYSTEM 

•  LANDING  GEAR 

•  ENGINE  FAILURE 

•  STORE  JETTISON 

•  LANDING  GEAR  RETRACTION 


LTV  Aarospaca  and  Dafansa  Company 

Aircraft  Division  PROGRAM  OPTIONS 


AUTOMATIC  WIND  OVER  DECK  SOLUTION 


TWO  CONSTRAINTS: 

•  MAXIMUM  SINK 

•  MAXIMUM  ANGLE  OF  ATTACK 

OR 

MAXIMUM  PITCH  RATE 
OR 

MINIMUM  LONGITUDINAL  ACCELERATION 


TWO  VALUES  DETERMINED: 

•  WIND  OVER  DECK 

a  STICK  DISPLACEMENT  {OR  TAIL  DEFLECTION) 
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LTV  Atrospae*  and  Dafansa  Company 
Aircraft  Olvialon 


PROGRAM  OPTIONS 


SOLUTION  TERMINATION 


•  POSITIVE  TMAX 

-  TIME  HISTORIES  STOP  AT  THE  SPECIFIED  TMAX 


•  NEGATIVE  TMAX 

-  TIME  HISTORIES  STOP  WHEN  A  POSITIVE  RATE  OF  CLIMB  HAS 
BEEN  ACHIEVED  AtjQ  ANGLE  OF  ATTACK  HAS  PEAKED. 

-  IF  A  POSITIVE  RATE  OF  CLIMB  IS  NOT  ACHIEVED  QB  ANGLE  OF 
ATTACK  IS  CONTINUOUSLY  INCREASING,  THE  TIME  HISTORY  WILL 
STOP  AT  THE  ABSOLUTE  VALUE  OF  TMAX. 


LTV  Aaroapaca  and  Dafansa  Company 

Aircraft  Division  PROGRAM  OVERVIEW 


AERODYNAMIC  DATA  INCLUDES 
PROPULSION-INDUCED  EFFECTS 

COEFFICIENTS  ARE  FUNCTIONS  OF: 

•  ANGLE  OF  ATTACK  OR  LIFT  COEFFICIENT 

•  NOZZLE  DEFLECTION 

•  THRUST  COEFFICIENT  OR  NPR 

•  FLAP  DEFLECTION 

a  TRIM  SURFACE  DEFLECTION 


RMCOOr 
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LTV  Aerospace  and  Defense  Company 
Aircraft  Division 


PROGRAM  OVERVIEW 


TYPICAL  AERODYNAMIC  DATA 
FOR  A  THRUST  VECTORING  CONFIGURATION 


AERODYNAMIC  COEFFXilENTS  HAVE 
THE  DIRECT  PROPULSION  EFFECTS 
REMOVED  AND  ARE  FUNCTIONS  OF: 

•  ANGLE  OF  ATTACK  OR  UFT 
COEFFICIENT 

•  THRUST  COEFFICIENT  OR 
NPR 

•  NOZZLE  DEFLECTION 

•  FLAP  DEFLECTION 

•  TRIM  SURFACE  DEFLECTION 


LTV  Aerospace  and  Defense  Company 

Aircraft  Division  PROGRAM  OVERVIEW 


FORCES  ACTING  ON  THE  AIRCRAFT 


in 
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AutoCAT  User's  Guide 


,  1 - 1 - 1 - 1 - 1 - 1 

4.  S.  e.  7.  8.  9.  10. 

TIMB  - 


Of  ATTACK  AS  A  FONCTJOH  OF  TIME 
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AIRBREATHING  BOOSTER  PERFORMANCE  OPTIMIZATION 
USING  MICROCOMPUTERS 

Ron  Oglevie 

Innne  Aerospace  Systems  Co. 

2001  Calle  Candela, 

Fullerton,  CA  92633 
(714)  526-6642 

AIAA  Astrodynamics  Conference 
Workshop  on  Trajectory  Optimization 

10  August,  1992 
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ft^lereeewipMtet  >eeed 

OVERVIEW 

Irvine 

OpMita, 

Aerospace 

Siwylilliii  •! 

Systems 

Tr^MMftee 

Company 

•  MICROCOMPUTER-BASED  OPTIMIZING  SIMULATION  OF 
TRAJECTORIES  (MOST)’ 

•  MOTIVATION  -  Fill  void  in  preliminary  design  tools 

•  Easy  to  use  fast  running  modes 

•  TPBV  solution  for  truth  model 

•  OTIS  PROGRAM  OPERATION  ON  PC 

•  LOW-THRUST  TRAJECTORY  OPTIMIZATION  PROGRAM  (MICROTOP) 


'  Work  performed  under  Air  Force  Contract  F33615-91-C-2100. 


Mlieeeiwewter  d 

MOST  TECHNICAL  OBJECTIVES 

Irvine 

Aerospace 

BEING  MET 

Systems 

Company 

•  SUITABLE  FOR  RAPID  PRELIMINARY  DESIGN 

•  MICROCOMPUTER  OPERATION  -  Run  time  less  than  5  mins,  on  PC  AT 

•  AIRBREATHING  &  ROCKET  PROPULSION  VIA  TABLES  AND 
EQUATIONS  -  Including  realistic  flight  constraints 

•  EARTH-TO-ORBIT  (ETO)  FLIGHT 

•  PLANAR  FLIGHT  -  Simple  rotating  earth  model  facilitates  3-D  type 
results  with  minimal  complexity 

e  EASE*OF-USE  -  Easy  input,  good  graphics,  &  robust  convergence 
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OpMwWl 


HYBRID  APPROACH  OFFERS 

Irvine 

Aerospace 

SPEED  AND  PRECISION 

Systems 

Company 

altitude  (FT)  -10^ 

62.8  '25.6  188.4  251.2  314.0 


1.  ^ _ _ 

2  STAGES  TO  ORBIT 

Irvine 

OpeuMMi 

Technologv 

SlWMlMlWI  Ot 

OPTIMIZATION  GOALS  ACHIEVED  ON  PC 

Group,  Inc. 

GRAPHICS  ENHANCE  UNDERSTANDING 
OF  PERFORMANCE  OPTIMIZATION 


Irvine 

Technology 
Group.  Inc. 


CONTOURS  OF  CONSTANT  MANEUVERING  EFFICIENCY 
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CONCLUSIONS  - 

IrvkM 

Aerospace 

Sbn^artanel 

ETO  PERFORMANCE  OPTIMIZATION 

Systems 

ACHIEVED  ON  PERSONAL  COMPUTER 

Company 

•  MOST  •  LOW  COST,  RAPID  RESPONSE  TOOL  FOR  PRELIM.  DESIGN 
SUCCESSFULLY  ACHIEVED 


•  ADVANTAGES  OF  PC  DEMONSTRATED  -  Low  Cost,  portabUHy,  and  good 
graphics  and  support  software  (LOTUS,  Harvard  Graphics,  Freelance,  etc.) 

e  USER-FRIENDLY  ENVIRONMENT  FOR  PREPARATION  OF  INPUT  FILES  &  OUTPUT 
DATA  •  Facilitates  OTIS  input  file  preparation 

•  MOST  FAST  RUNNING  MOOES  DEMONSTRATED  -  Good  agreement  with  OTIS 
results.  Early  engineering  model  delivered 

•  PRELIMINARY  RESULTS  FROM  2-D  NLP/COLLOCATION  ALGORITHMS  (MINI- 
OTIS)  ARE  ENCOURAGING 

•  FAST  RUNNING  MODES  FACILITATE  NEW  APPLICATION  -  Trajectory  optimizer 
simple  and  fast  enough  to  imbed  in  vehicle  design  optimization  code 

•  OTIS  HOSTED  ON  PC  •  ETO  flight  acNevable  with  large  RAM  ( ~40KBYTES) 
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An  Algorithm  for  Trajectory  Optimization  on  a 
Distributed-Memory  Parallel  Processor 

Mark  L.  Psiaki  and  Kihong  Park 
Sibley  School  of  Mechanical  and  Aerospace  Engineering 
Cornell  University 


Acknowledgement: 

Work  supported  by  NASA/LaRC 
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Continuous-Time  Problem  to  be  Solved 


find:  i<(t)  and  jc(t)  for  to  ^  t  <  tf 

tf 

to  minimize:  ^  “  J  dt  +  V[jc(tf)] 

to 

subject  to:  jc(to)  given 

i  =yij:(t).tt(t),t] 
ae[j:(t)>ii(t),t]  =  0 
aiWt),M(t),t]  <  0 
«ef[4tf)]  =  0 

flifWtf)]  ^  0 


Approach 

•  Use  zero-order-hold  control  parameterization 

•  Model  as  a  multi-stage  parameter  optimization  problem 

•  Retain  state  variables  and  dynamic  constraints  explicitly 

•  Solve  using  a  nonlinear  programming  algorithm  that ... 

...  has  fast  local  and  robust  global  convergence 

...  allows  infeasible  intermediate  results 

...  parallelizes  function,  gradient,  etc.  evaluations  at  different 
time  steps 

...  exploits  dynamic  structure  and  parallelism  to  get  search 
directions 
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Multi-Stage  Nonlinear  Programming 

Problem 

■  T  T  T  T  T  T  I- 
^  =  ^0*  *  1»  ^1*  ^  ^  N-l» 

N-1 

to  minimize:  J  =  X  Lk(*k»“k)  +  ^[xnI 

k=0 

subject  to:  Xq  given 

*k+i  =  4(*k»Uk)  for  k  =  0  ...  N"1 
aeit(X|j,U|f)  =  0  for  k  =  0  ...  N- 1 
ai,^(X|(,Ui()  <  0  fork  =  0...  N-1 
aeN(*N) =  0 

aiN(XN)  <  0 


A  Static/Dense  Nonlinear  Programming 

Problem 

find:  X 

to  minimize:  J  (x) 

subject  to:  Ce(x)  =  0 

Ci(x)  <  0 
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Time  to  Solve  QP  on  the  INTEL  iPSC/2  (sec.) 


Status  of  Project 

•  Parallel  search  direction  algorithm: 

FORTRAN  version  tested  on  32-node  INTEL  iPSC/2 

•  General  static  NP  algorithm: 

FORTRAN  version  tested  on  1  node  of  INTEL  iPSC/2 
Compared  to  NPSOL  version  4.02  on  static  problems 

•  Full  parallel  trajectory  optimization  algorithm: 

A  "next  generation"  of  the  NP  algorithm  that  exploits 
parallelism  and  dynamic  problem  structure 

FORTRAN  components  currently  being  tested  on  32-node 
INTEL  iPSC/860 
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Plans 

(the  LORD  willing) 

•  Finish  component  and  full  algorithm  testing  (Present-Oct.  ’92?) 

•  Model  and  solve  guidance  problems  for  NASP  and  generic 
hypersonic  vehicles  (Oct.  '92  -  Dec.  '93) 

•  Compare  to  existing  codes  (199?) 

•  Evaluate  suitability  for  real-time  guidance  updates  (199?) 

•  Make  code  user-friendly  and  disseminate  (199?) 


DIstrlbytlon  of  Problem  Sieges  on 
Pereliel  ProceSSoro 


24  Stoae  Probigm  on  8  Proccsgori 


0.1.2 


101 


Divide-and-Conquer  Trajectory 
Optimization 

Iteration  I  x  is  fixed  at 


Iteration  2  x  is  fixed  at 


X  is  flxed  at 
this  time 


Iteration  4 


t 


Test  Problem  2  pC'P  Si  ATX  C  A^P 
find:  xi,  X2 

to  minimize:  J  =  -  X2 

subject  to;  (xj  -  l)^  +  X2  lOCXX)  (X|  X2  -  1)^  *  .0625  <  0 


Test  Results 

•  Augmented  Lag.:  77  iterations 

(52  feasibility  and  25  optimality) 

•  NPSOL  4.02;  feuledi 
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Test  Problem  S  FcO.  STATTC  MP 
find:  Oq,  AVt.  Av,  AV2 

to  minimiro:  J  s  Oq 

subject  to:  Newton's  laws  for  a  spherical  Eanfa 

Fixed  fuel  specific  impulse 

ri3)-er 

Vc«.-ev  i  V,  i 
-  6^  ^  y,  ^  +  e, 

28o  -  ej  i  if  i  28®  q 
°>a^xy  ^  ®f 


Test  Results 

Augmented  Lag.:  iterttfVeni 

»/Ma  '  . 

n/i^ 


^ZMXr  ■  /oo  ^ 

•  NPSOL  4.02:  9  iterations 


Aero-Assisted  Orbital  Maneuvering 

Example 

(taken  from  MielCt  1989  ACC) 

Problem:  Minimize  Fuel  for  GEO  to  LEO  transfer  with  +28^  inclination 
change 

x  =  [v,  Y,  \|/,  r,  <|),  e]^ 
u  =  [Vc,  Yc,  VcF  or  [Cl.  o,  t]’’ 

Constraint:  Heating  rate  ^  1(X)  watts/cm^ 

LOR-like  problem  derivation:  Linear-quadraticize  about  guessed  solution 
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GEO 


1 


